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1. Introduction 

Confinement of quarks, at feast on the fattice seems to be due to the formation of a fiux tube be- 
tween a quark and an anti-quark in the QCD vacuum [|]]. It has been conjectured that the properties 
of tliis flux tube can be described by an effective hadronic string 

On the lattice, this flux tube can be observed by measuring the potential between a static quark 
and an anti-quark. One of the characteristics of the string like behaviour of the flux tube is the 
presence of a long distance 1 /r term in the qq potential in all dimensions and with a universal 
coefficient. It is known as the Luscher term [^, ^]. 

The Liischer term has been looked at in lattice simulations since the eighties [||, ^. However 
in recent times, increase in computing power and improvement in algorithms have allowed really 
precise measurements of that term and also the next subleading r^^ term. See [0-[T3|] for example. 

Another aspect of the string behaviour is the excitation spectrum of the flux tube. A recent 



review on this topic can be found in [14| 



In this article we present results for the continuum limit of our simulations of the Polyakov 
loop correlators for = 3 SU (2) Yang-Mills theory and compare the resulting static potential with 
both perturbation theory and string model predictions. This allows us to narrow down bounds on 
the distance beyond which we can say the flux tube indeed shows a string like behaviour. 



2. Simulation parameters 

We carried out simulations of three dimensional SU (2) lattice gauge theory on lattices at four 
different lattice spacings. On these lattices the scale was set by the Sommer parameter ro = 0.5 fm, 
implicitly defined by r^firo) = 1.65, where f{r) is the force between the static quark and the 
anti-quark. Our measured values were ro/a = (3.9536(3), 6.2875(10), 8.6602(8), 10.916(3)) at 
j3 = (5.0,7.5, 10.0, 12.5) respectively. We also observed that ar^ (a : string tension) is a constant 
(~ 1.522) to a very good approximation as expected in the continuum limit. With this scale our 
coarsest lattice had a spacing of slightly below 0.13 fm and our finest lattice spacing is about 
0.045 fm. We used symmetric cubic lattices and the Wilson gauge action. On all these lattices, 
we computed Polyakov loop correlators {P*{x)P{y)) for various spatial separations r = y — x, with 
separations being taken along the axes only. 

To reliably extract signals of these observables which are exponentially decreasing functions of 
r and T (the temporal extent of the lattice), we used the Luscher- Weisz exponential error reduction 
algorithm [|8]|. For further details of our run parameters, see [p3[]. 



3. Results 

From the {P*P) correlator one can extract the static quark-antiquark potential V{r) by V{r) = 
— ^ln(P*/'(r)). We prefer to look at the first and the second derivative of this potential for the 
force between the quark and the antiquark and information about sub-leading terms. To facilitate 
our comparison with string models we actually compute a scaled second derivative which we call 
c(r). This quantity is expected to become the Liischer term (= —^^-^^) asymptotically. 
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Figure 1: r^f{r) vs r/ro. The horizontal line locates the Sommer scale. 



On the lattice these are given by 

F{f)=V{r)-V{r-\) & c{r) = j[V {r+ I) + V{r - \) -2V{r)] (3.1) 

where r = r + | + &{a^) and r = r + G{c^^ are defined as in [Q] to reduce lattice artifacts. 

For large distances, we are going to concentrate on the potential due to the NG string, the 

so called Arvis potential [ p^ given by VArvis = err ^1 — ^jj^jr^ ■ We will define the L.O. and 
N.L.O. approximations by retaining 1/r and 1/r^ terms in the potential respectively. We will 
compare our lattice data on force and c(r) with predictions from leading order, NLO and from the 
full Arvis potential. 

At short distances we compare with the perturbative potential obtained by Schroder [ [F^ ] as 
Vpert(r) = sp,nr+ ^ lngV + . . . with ^rt = For SU{2) Cp = 3/4, Ca = 2. The perturba- 

tive force and Cpert(?-) can be computed by /pert(?") = "^^^^ and Cpert(?-) = ^-^p^. 

In Fig. [l] we plot r^f{r) versus r/vQ. The horizontal line is r^f{r) = 1.65 and defines the 
Sommer scale tq. The N.L.O. curve lies in between the L.O. and Arvis curves and has been omitted 
for the sake of clarity. The data starts departing from the one-loop perturbative curve around 0.22 
ro or 1.8 GeV and joins onto the string curves around 1.5 ro or 260 MeV. The scaling exhibited by 
the data is very good with all the four different beta values falling on the same curve. Beyond 1.5ro 
it is virtually impossible to distinguish the different theoretical curves as they are all dominated by 
the string tension. The force data in fact gives the impression of the string description being good 
even at distances as small as tq. 

c(r) does not contain the string tension, and has a universal value in the L.O. It is therefore 
more sensitive to the sub-leading behaviour of the flux tube. In Fig. ^ we plot c(f) in units of r/ro. 
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We look at a wide range of r starting from where the data almost touches the perturbative curves 
going all the way to the region where the string predictions hold. 

The data almost lie on top of each other exhibiting nice scaling behaviour as one goes to larger 
values of r. The j8 = 12.5 and j8 = 10 data come together already in the range 0.5 and 1.25 tq. The 
j3 = 7.5 set joins onto this at around 1.5ro and even the j8 = 5 data joins up at around 2.25ro. This 
points to the possibility that the continuum limit of the scale where the flux tube is well described 
by the Arvis curve can be obtained even on relatively coarse lattices. 

We estimate a temporal extent correction factor of about 0.1% for c(r) for our largest r values. 
This is about an order of magnitude lower than our statistical errors at such r values. Corrections 
due to finite spatial extents are of similar magnitude. 

Ind = 3, IR divergences prevent computation of the perturbative potential beyond one loop'. 
We tried to obtain these terms from our data. Assuming the perturbative potential to be of the form 

Vpert(r) = ^ln/r+Z^^r+AAlngV + BgV + ... (3.2) 
271 6471 

we obtain A = 0.013162(3) and B = 0.001089(1). Also we estimate the range of validity of first 
order perturbation theory to be about 0.1 fermi (consistent with our estimate from the force data). 

From our data we find that the string tension constitutes 95% of the force at around 1.3ro, 98% 
at around 2.1ro and 99% at around 2.9ro. The relative difference between the Arvis and the leading 
order force, which gives us an idea about the importance of the subleading behaviour, is about 2% 
at 1.02ro, 1% at 1.2ro and goes down to 0.1% at about 1.9ro. 

'NDH wishes to thank G. 't Hooft for an illuminating discussion on this. 
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r 


r 


P=5 


P 


= 7.5 


j3 = 10 


p = 12.5 


3 


2.808 


-0.10274(3) 






-0.06172(1) 


-0.05095(1) 


4 


3.838 


-0.11886(8) 






-0.07737 (3) 


-0.06530(3) 


5 


4.876 


-0.1265 (2) 


-0 


1052(1) 


-0.08937 (6) 


-0.07714(5) 


6 


5.903 


-0.1299 (4) 


-0 


1128 (2) 


-0.0985 (1) 


-0.0868 (1) 


7 


6.920 


-0.1322 (9) 


-0 


1183(3) 


-0.1060 (5) 


-0.0946 (2) 


8 


7.932 


-0.1342(10) 


-0 


1222 (3) 


-0.1118 (7) 


-0.1011 (2) 


9 


8.941 


-0.1333 (6) 


-0 


1251 (6) 


-0.1156 (2) 


-0.1064 (2) 


10 


9.948 


-0.1340(10) 


-0 


1273 (8) 


-0.1190 (3) 


-0.1108 (3) 


11 


10.95 


-0.1375 (34) 


-0 


1288 (11) 


-0.1217 (4) 


-0.1145 (4) 


12 


11.96 




-0 


1296 (7) 


-0.1241 (5) 


-0.1174 (5) 


13 


12.96 




-0 


1308 (9) 


-0.1258 (7) 


-0.1201 (6) 


14 


13.96 




-0 


1323 (9) 


-0.1270 (6) 




15 


14.97 




-0 


1332(13) 


-0.1282 (9) 




16 


15.97 




-0 


1345 (21) 


-0.1288 (11) 




17 


16.97 








-0.1303 (15) 




18 


17.97 








-0.1308 (23) 





Table 1: c(r) 



The type of the string is determined by the sub-leading behaviour of the flux tube as at leading 
order, a variety of theories yield the universal Liischer term [18, 19 1. It is clear from the data that 
the approach to the Liischer term is from below, consistent with effective bosonic string model 
predictions while at short distances the data matches perturbation theory. Beyond a distance of 
about 2.75ro, the data seems to be well described by the Arvis curve. 

An interesting question is what happens to this distance for SU {N) as N increases. For SU (3) 
on coarsest lattices it is about 1.8ro. Similar indications are also there for SU (5) [|I2[]. However 
this scale shifts towards larger r as one approaches the continuum limit. 

Beyond N.L.O., effective string theories parametrise c(r) to be of the type a(ro/r)'* + j8(ro/r)^. 
Our data gives a = 0.209(9) and jS = —0.235(24), which is quite different from the Arvis values. 
It is therefore interesting to know the effective string predictions for these coefficients. 

At intermediate distances over a wide region of r varying from 0.5ro to 2.8ro, the data is very 
well described by a formula of the type 

Existing SU (3) data in 3-d also admit a similar description. The curves are shown in Fig. |3[ It 
is not clear if there is any theoretical basis for such a description, but at least they provide accurate 
interpolation formulae. 



4. Conclusions 

In this article we have looked at the continuum limit behaviour of the SU (2) flux tube at 
intermediate distances by measuring the static qq potential. Starting from a distance of about 0. 1 
fermi where the potential starts breaking away from 1-loop perturbation theory, we go to distances 
of about 1 .4 fermi where the data begins to be well described by the Arvis potential. 
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Our data on c(f) seems to approach the Liischer term from below as expected in bosonic string 
models. At distances below 0.15 fermi the data joins onto the perturbative values. An empirical 
formula describes the data well in the intermediate region. 
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